Nonlocal PNJL model beyond mean field and the QCD phase transition 
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A nonlocal chiral quark model is consistently extended beyond mean field using a strict 1/N C 
expansion scheme. The parameters of the nonlocal model are refitted so that the physical values of 
the pion mass and the weak pion decay constant are obtained. The size of the 1/N C correction to 
the quark condensate is carefully studied and compared with the usual local Nambu-Jona-Lasinio 
model. It is found that even the sign of the corrections can be different. This can be attributed to the 
mesonic cut-off of the local model. The model is also applied to finite temperature. We find that the 
1 /N c corrections dominate the melting of the chiral condensate at low temperatures, T < 100 MeV, 
in agreement with chiral perturbation theory. On the other hand, the relative importance of the 1/N C 
corrections in the cross-over regime depends on the parameter To of the Polyakov loop potential. 
For To = 270 MeV, corresponding to a fit of lattice data for pure gluodynamics, the correction 
terms are large and lead to a lowering of the chiral phase transition temperature in comparison 
with the mean-field result. Near the phase transition the 1/N C expansion breaks down and a non- 
perturbative scheme is needed to include mesonic correlations in that regime. Lowering To leads to 
a more rapid cross-over even at the mean- field level and the unstable region for the 1/N C corrections 
shrinks. For To J$ 220 MeV the temperatures of deconfinement and chiral restoration are practically 
synchronized. 
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I. INTRODUCTION 

A quantum field theoretical description of strong inter- 
actions in the nonperturbative regime is one of the most 
interesting and challenging problems of present-day theo- 
retical physics. Quantum chromodynamics is well known 
only at the perturbative level whereas the low-energy re- 
gion and the most interesting "hadronic" phase in the 
QCD phase diagram is in the nonperturbative regime. 
The only nonperturbative ah initio calculations are per- 
formed in lattice QCD, but their range of applicability 
is still limited. To gain some analytical insights to non- 
perturbative QCD, continuum approaches, even using ef- 
fective models, are legitimate tools. They may provide 
a theoretical interpretation of results from lattice QCD 
and allow their extrapolation to otherwise inaccessible 
domains. 

One of the successful models for a description of chiral 
quark dynamics and the phase diagram is the Nambu- 
Jona-Lasinio (NJL) model [l[ applied to quarks 
This model provides a mechanism for spontaneous chi- 
ral symmetry breaking and the formation of a quark 
condensate. The low-lying hadron spectrum, low-energy 
dynamics, the main strong and electromagnetic decays, 
hadron-hadron scattering and the internal characteris- 
tics of mesons have a reasonable explanation within this 



model. A generalization of the NJL model has been pro- 
posed which includes the coupling of the chiral quark sec- 
tor to the Polyakov loop, being an order parameter of the 
deconfinement transition (PNJL model fa-tl5|). The tem- 
perature dependent parameters of the effective Polyakov 
loop potential arc fitted by using lattice QCD data on 
pure gluodynamics. 

Nonlocal generalizations of the PNJL model provide an 
approach to the 4-momcntum dependence of the quark 
mass function and wave function renormalization of the 
quark propagator [H|-[l3], thus allowing to implement 
detailed nonperturbative information about low-energy 
QCD dynamics accessible, e.g., in ab-initio LQCD sim- 
ulations fl8l - l2C| . Recently, a justification of nonlocal 
PNJL models as effective low-energy limit of QCD has 
been given 17J ) based in part on methods of the Wilso- 
nian renormalization group (2lj . 

Usually, NJL and PNJL models are formulated at the 
mean-field level. However, there arc physical problems 
where the mean-field formulation is not sufficient. Large 
corrections to the mean-field behavior can be expected, 
e.g., in the description of broad resonances from their 
coupling to intermediate meson states 1 and for the equa- 
tion of state of the hadronic phase where quark and 
gluon degrees for freedom are "frozen" in condensates 
and hadronic bound states are carrying all dynamics and 
thus represent the physical degrees of freedom. 
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The cr-meson with a decay width of the order of its mass is the 
most striking example of such state. See, e.g., "Note on scalar 
mesons" in 12211. 
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There are different schemes to go beyond the mean- 
field level [23|-[33j]. One of the most promising ones is 
based on a strict expansion of the inverse number of 
quark colors, 1/N C , which is a natural expansion param- 
eter for gauge theories 34| . The local N JL model is non- 
renormalizable and therefore it is necessary to introduce 
an additional cut-off parameter when going beyond the 
mean-field level. This problem is absent in nonlocal ver- 
sions of the NJL model where the nonlocality leads to an 
effective regularization which renders the quark (multi- 
loop diagrams convergent. 

In the present paper the SU(2) x SU(2) nonlocal PNJL 
model is investigated beyond mean field within a strict 
l/N c expansion scheme at finite temperature and zero 
chemical potential. This study aims to quantify the effect 
of mesonic excitations on the chiral restoration tempera- 
ture and to demonstrate that the behavior of the chiral 
condensate at low-temperatures is in accordance with the 
exact results of chiral perturbation theory. We also in- 
vestigate the dependence of the results on the parameter 
Tq of the Polyakov-loop potential. 



II. NONLOCAL MODEL IN VACUUM 

We begin with the discussion of our nonlocal model in 
vacuum. The main goal is to fix the model parameters by 
calculating meson properties at next-to-lcading order in 
l/N c . It is well-known that in the local PNJL model at 
T = the gluon sector decouples from the quark sector, 
so that the latter is reduced to the standard NJL model. 
The same is true in the nonlocal model. In this section 
we can therefore restrict ourselves to the quark sector, 
while the Polyakov-loop dynamic will be introduced in 
Sec. 



A. Mean field 

The quark sector of the nonlocal chiral quark model is 
described by the Lagrangian 

C q = q(x)(i$-m c )q(x) + ^[J 2 (x) + J 2 (x)], (1) 

where m c is the current quark mass. The nonlocal quark 
currents are 

J M (x) = J d i x 1 d i x 2 f(xi)f(x 2 ) x 

xq(x- Xi)T M q(x + x 2 ), (2) 

where T a — 1, IV = ij 5 r a with a = 1, 2,3, and f(x) is 
a form factor. The latter is defined by its Fourier trans- 
form in Euclidean space 2 , which we take to be Gaussian, 



f 2 (p%) = exp(— p 2 E /K 2 ). This scheme introducing non- 
local currents which we utilize throughout this work has 
been denoted as "instanton liquid model" as opposed to 
the "one-gluon exchange model" scheme, see Ref. [35| for 
details. 

After linearization of the four-fermion vertices by intro- 
ducing auxiliary scalar (a) and pseudoscalar (7r a ) meson 
fields the quark sector is described by the Lagrangian 3 



C 



q{x)(i@ - m c )q(x) 



n 2 a (x)+a 2 (x) 



2G 



+J a (x)a{x) + ir a (x)J%{x) 



(3) 



To proceed, we single out the nonzero mean-field value of 
the scalar field by the decomposition a = a + ctmf so that 
7r a and a denote only the fluctuating parts of the fields 
((ir a ) = (a) = 0) describing mesonic correlations. The 
scalar mean field gives a dynamical contribution to the 
quark mass, i.e., the dressed quark propagator becomes 

-l 



S(p) = {j>-m(p 2 )y 
with the quark self-energy 



0- 



(4) 



iGT a {jp,p) 

m d f 2 (p 2 ) 



d 4 k 



Tr [T a (k,k)S(k)] 



(5) 



Here 
color- 



the symbol 
flavor- and 



Tr stands for the trace c 
Dirac-indices, and F M (qi,q 2 ) 



^Mf{q 2 )f(q 2 ) is the nonlocal generalization of the 
meson-quark-antiquark vertex function with the quark 
and antiquark momenta q\ and q 2l respectively. In the 
following, we will often use a shorthand subscript no- 
tation for the momentum dependence of functions, e.g., 
f k = f(k 2 ) and T-( qi ,q 2 )=r^ q2 . 

The amplitude = -<tmf is an order parameter for 
dynamical chiral symmetry breaking. The chiral conden- 
sate per flavor, 

(m) MF = 4 f^z^[s np m (6) 



(2nY 



is obtained from the non-perturbative part of the quark 
propagator, S np (p) = S(p) — S c (p), i.e., after subtracting 
the perturbativc part S c (p) = (J> — m c ) . 

Mesons are described as bound state solutions of the 
quark-antiquark Bethc-Salpeter equation. The meson 
propagators are given by 

m 



D M = 
v 



where M 



7T, a and Uf 



-G- 1 +IJM' 

IIy F (p 2 ) are the mean field 



polarization functions defined by 
d 4 k 



n 



M 



(2ny 



Tr 



S k _T 



M 



_dk + L k+,k- 



(8) 



2 Unless stated otherwise, the expressions in this paper are given 
in Minkowski space. The transformation to Euclidean space is 
trivial for energies below the (pseudo)threshold (see Sect. Ill Ct . 



3 Note, that after bosonization the possible exchange (Fock) terms 
between quark currents are eliminated. 
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where k± = k ±p/2. 

The meson masses are the poles of the propagators at 



P~ 



(M^ F ) 2 obtained by solving 



G~ 



ttMF 



0. 



(9) 



Let us consider the pion case. In the vicinity of the pole 
the pion propagator can be expanded as 



v 



(5. MF ) 



p 2 _ (MMF)l 



+ regular terms , 



(10) 



where g% is the pion-quark-antiquark coupling constant 



I MF\ 
(&r ) 



dp 2 



(11) 



p 2 = ( M MF)2 



To calculate the weak pion decay constant, the pio- 
nic solutions of the Bethe-Salpeter equation have to be 
coupled to an external weak current. To that end the 
Lagrangian ([3]) must be modified so that it becomes in- 
variant under local vector and axial- vector gauge trans- 
formations. In nonlocal models this is complicated by the 
fact that not only the kinetic part part but also the inter- 
action is not gauge invariant by itself, so that in addition 
to the usual covariant derivative one has to take into ac- 
count the coupling of the external fields to the nonlocal 
quark vertices. This can be done by derealization of the 
quark fields 



q{y) Q{x,y) = E(x,y)q(y) 



(12) 



where 



E(x, y) = Vexp { i \ dz»[V a Jz) + A a Jz) l5 ]T a } (13) 



is the Schwinger phase factor, involving the external 
vector and axial- vector gauge fields V° and *4^, and 
T a = r a /2. For the kinetic part it can be shown that 
this replacement is equivalent to minimal substitution, 
hence leading to the standard electroweak vertices. In 
the interaction Lagrangian, on the other hand, the non- 
local quark currents Eq. @ arc replaced by 



Jm{x) 



d 4 x 1 d 4 x 2 f(x 1 )f(x 2 ) x 

xQ(x - xi,x)r M Q(x,x + x 2 ), (14) 



giving rise to additional vertex contributions. One of 
them is related to the coupling of J a to the scalar mean- 
field ctmf and contributes, together with the bare vertex 
extracted from the kinetic part, to the vertex function 
depicted on the left of Fig. [TJ The other contributions 
originate from Jm coupled to the fluctuating meson fields 
and lead to the vertex functions shown on the right of 

Fig.ru 

Although the above procedure is sufficient to ensure 
gauge invariance, in general it does not unambiguously fix 





FIG. 1: Vertices of weak currents (wavy lines) coupled to a 
quark (left) and to the non-local quark- meson vertices (right). 





1 2 
FIG. 2: Weak pion decay at mean field. 



the vertex structure. To this end, it is necessary to define 
rules for the evaluation of the line integral in Eq. (|13p . 
This can be done by specifying the integration path, e.g., 
as a straight line [35|, 1371 l38l . |4fJ , or by making use of the 
path independent definition of the derivative of the line 
integral [HI, H^|. However, for the pion decay constant 
we only need the longitudinal projection of the vertices, 
which are related to axial Ward-Takahashi identities and 
therefore do not depend on the integration path. The 
results read 



k— ,&4 



1 



S^ l5 T a + l5 T a S^- 



+ 



+2(m c + m d f k _f k+ ) l5 T a j , 

(fk+ - fk + - P ) fk-lsT a T M + 

(fk- — fk-+p) /fe + r M 757 ia 



(15) 



Accordingly, the pion decay constant at the mean-field 
level contains two pieces, f™ F (p 2 ) = f™J (p 2 ) + f™J (p 2 ), 
which are displayed in Fig. fJl Evaluating these diagrams 
one finds 



/5V) 
/. M 2 V) 





f d 4 k 


Tr 


1 k + ,k 




f d 4 k 
(2tt)4 


Tr 


1 k,k 



k + ,k-Sk_ r ' k _ yk+ dk + 



5,L 



Sk 4 



(16) 



With the definition of the pion mass and weak decay 
constant two of the three parameters (m c , A, GA 2 ) of 
the nonlocal model can be fixed with the vacuum values 
of these observables. 



B. 1/N C corrections beyond mean field 

As usual in the systematic \/N c expansion, the four- 
quark coupling constant G is considered to be of the order 
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FIG. 3: 1/N C corrections to the quark selfenergy. 



1/N C . In this case the N c behavior of pion properties cal- 
culated in the model coincides with leading-order QCD 
calculations (/ w ~ y/Nc)- As a result any meson propa- 
gator line in diagrams has a 1/N C suppression factor. 

The next-to-leading 1/N C correction E N ° to the quark 
selfenergy corresponds to the diagrams displayed in 
Fig. [3] [27H29I | . From these diagrams one obtains 

v Nc _ n ( 2 \ "* „■ f ^ TrvM^M o r M 1 

^-p — W P Z-' J (2?r)4 L 1 1 pv-i a P~ ll p-i,p\ > 



M 



where 



c=djc, c= E 1 



M = CT.7T 



dH 



-pjMpcrMM 



(17) 



and ^qfq 2 is the quark triangle diagram for the three- 
meson vertex. For external mesons a, j3 and 7 the quark 
triangle has the following form 



91,92 



j£k_ 



Tr x 



(18) 



1 k+qi ,k+q 2 



When qi = 52 = Z, the three- meson vertices in Euclidean 
space are 



-pcrMM 

1 U 



-4N r N 



d%k fkfk+i 



x [2* • (fc + l)m k+l ± m fe ((fc + If 



(19) 

a, the lower 
m 2 is the denominator of 



where the upper sign corresponds to M 
sign to M = 7r and Dk = k 2 
the quark propagator. 

From the 1/N C corrections to the quark selfenergy it 
follows for the quark propagator 



qMF+N c 



S n 



Sp^p c Sp 



higher orders, (20) 



which, in turn, gives rise to 1 /N c corrections to the quark 
condensate (cf. Eq. ([6])). According to the two selfenergy 
corrections shown in Fig. [3j there are two contributions, 
(qq) Nc = (qq) Nc ' a + (qq) N "' b , which are given by 



(qq)^ a = -Clir 



(qq) N -» = ( b 



(21) 



where lip' can be obtained from Eq. ([8]) for IEJ by sub- 
stituting the form factors in the numerator fl f?_ — > 



fk + fk_ an( i C b can be obtained from £ by substituting 
the form factors in the numerator of the function rf f lM , 

Eq. m, fin ft+i ■ 

The full pion propagator consists of a mean-field part 
plus 1 /N c corrections (details are given in the appendix) 

nr(p 2 ) = nf( P 2 ) + n^( P 2 ) (22) 

and the pion mass can be found by solving the equation 

- G- 1 + nl ull {(M™ 1 ) 2 ) = 0. (23) 

The quark-meson coupling constant should be divided 
into mean- field part and a 1/N C correction. This sep- 
aration is unique only for the pion in the chiral limit. 
At finite m c there are different possibilities to take into 
account higher order p 2 terms. In the present work we 
use 



go 
(9 



dp 2 
dnf F ( P 2 ) 



dp 2 



p2 = ( A fFull)2 



p 2 = ( M MF)2 



(24) 



There are two sources oil/N c corrections for the weak 
pion decay constant. One is the correction of the meson- 
quark coupling constant and can be obtained from the 
mean-field expression, Eq. (|16p . by the substitution 



g^ F — > g^ c = g Fu11 — g^ F ■ The other correction is due to 
new diagrams appearing at order 1/N C . We present de- 
tails of the calculation of this correction in the appendix. 



C. Vacuum results 

The 1/N C corrections to meson properties will affect 
the results for the quark condensate via the readjust- 
ment of the model parameters (A, m c , GA 2 ) which are 
to be chosen such that the physical values for the pion 
mass M v ± = 139.57 MeV and the weak pion decay con- 
stant / w = 92.42 MeV are obtained at T = 0, while 
the dimcnsionlcss coupling GA 2 is left as a free param- 
eter. Different parameterizations of the nonlocal model 
beyond mean field arc given in Tabic |TJ The correspond- 
ing behavior of the quark condensate as a function of the 
dimensionless coupling is presented in Fig. [4] 

One important check of the calculations is the proof of 
the Goldstone nature of the pion. In the case of exact 
chiral symmetry the pion should be masslcss. This check 
is shown in Fig. [5] where M Ful1 and M^ F are displayed 
as functions of the current quark mass for A and G taken 
from parameter set 4. 

The mean-field contributions to the pion mass and 
weak decay constant are also listed in Table Q] for the 
different parameter sets. The corresponding 1/N C cor- 
rections are shown in Fig. [6j For lower values of GA 2 
these corrections amount to —15 MeV for the pion mass 
and to 20 MeV for f v . For parameter set 4 the correc- 
tions are only about —2 MeV and 5 MeV, respectively. 
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TABLE I: Different parameterizations: model parameters (A,GA 2 ,m c ), the dynamical mass rrid, position of the first two poles 
and the corresponding (pseudo-) threshold of the mean-field quark propagator, mean-field values of the pion mass and the 
weak decay constant, and critical temperatures for chiral restoration and deconfinement for two values of the parameter To of 
Polyakov loop potential. The set with the highest threshold value is highlighted. The pseudo-threshold values are shown in 
italics. 



No. 
units 


GA 2 


A 

MeV 


m c 
MeV 


TCld 

MeV 


first two poles 
GeV 2 


(pseudo)threshold M™ F 
MeV MeV 


fMF 

MeV 


To = 


270 MeV 


To = 


208 MeV 


T c 


rpMF 

* c 

MeV 


T d 


T c 


c 

MeV 


T d 


f 


13.35 


1479.2 


2.82 


139.2 


-0.0205, -6.331 


286.7, 5032.3 


155.5 


71.5 


191 


205 


213 


164 


162 


162 


2 


14.89 


934.8 


5.58 


211.2 


-0.0529, -1.545 


459.8, 2486.0 


144.6 


83.4 


193 


206 


214 


166 


165 


166 


3 


17.06 


705.9 


8.64 


269.1 


-0.1262, -0.447 


710.6, 1337.2 


142.4 


87.1 


193 


207 


214 


168 


167 


168 


4 


17.64 


670.3 


9.38 


281.9 


-0.2291 ±i0 


957.3 


142.2 


87.6 


194 


207 


214 


169 


168 


169 


5 


19.72 


580.5 


11.78 


322.5 


-0.1082 ± i0.1763 


793.7 


141.7 


88.7 


197 


208 


213 


170 


169 


170 


6 


22.83 


500.8 


14.95 


373.8 


-0.0289 ±i0. 1832 


654.7 


141.4 


89.5 


200 


208 


213 


170 


170 


170 


7 


26.33 


445.3 


18.15 


424.0 


+0.0117 ±i0. 1692 


562.0 


141.2 


90.0 


202 


208 


212 


171 


171 


171 


8 


30.20 


404.4 


21.37 


473.4 


+0.0344 ±i0.1536 


496.0 


141.2 


90.3 


203 


208 


211 


172 


171 


172 



280 

260 
% 240 
co 220 
^ 200 
T 180 

160 

140 

14 16 18 20 22 24 26 28 30 

GA 2 

FIG. 4: (Color online) The quark condensate as a function 
of the dimensionless coupling GA 2 . The result of the full 
1/N C corrected approach (black solid line) is compared to the 
mean field contribution (red dashed line). The vertical line 
corresponds to the parameter set No. 4. 
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FIG. 5: (Color online) Squared pion mass normalized to the 
physical value versus normalized current quark mass in mean- 
field approximation (red dashed line) and the full 1/N C cor- 
rected approach (black solid line). The values of A and G 
correspond to parameter set No. 4. 



The region of applicability of the model is closely re- 
lated to the pole structure of the mean-field quark prop- 
agator and depends on the parameters. For large values 
of A, the propagator has two poles on the real p 2 axis, 
whereas when A is decreased these poles eventually merge 
and then go over into a pair of complex conjugate poles. 
(In both cases, there is an infinite number of additional 
complex conjugate poles, which, however, only play a 
minor role.) If the propagator has real poles, there is a 
real threshold in quark loop diagrams, i.e., for external 
momenta p > ^threshold the loops have imaginary parts. 
The complex conjugate poles, on the other hand, lead to 
a so-called pseudo-threshold. In the absence of real poles 
the quark loops are then purely real but have a cusp at 
the pscudothrcshold. Usually, the absence of poles on 
the real p 2 axis is considered as a possible criterion for 
quark confinement 41(. However, due to the cusp in the 



real part of the quark loop diagrams the applicability of 
the model above the pseudo-threshold is at least ques- 
tionable. 

The positions of the two lowest quark poles and the 
corresponding (pseudo-) thresholds of the polarization 
quark loop are listed in the table. The (pseudo-) thresh- 
olds are also shown in Fig. [7] There one can see that 
parameter set 4 corresponds to the maximal threshold 
value. For this reason we will use parameter set 4 in the 
finite T calculations in Sec. IIIII 



D. Sign of the 1/N C corrections 

Our results for the l/iV c -corrections to the quark con- 
densate and to the pion decay constant are quite surpris- 
ing. Naively, one would expect that taking into account 
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FIG. 6: (Color online) 1/N C corrections to the pion mass 
(black solid line) and the weak pion decay constant (red 
dashed line). 
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FIG. 7: (Color online) The (jg-(pseudo)threshold due to the 
lowest singularities of quark propagator. Below a critical 
value of GA 2 (thin solid line) there are real thresholds in 
quark loop diagrams (green dotted line for the first singu- 
larity and red dashed line for the second one). Above this 
critical value the lowest singularities of the quark propagator 
are complex-conjugated values and there is no imaginary part 
of the quark loops. However, quark loops contain a pseudo- 
threshold (thick solid line). The blue dash-dotted line indi- 
cates the cr-meson mass at the mean field level. 



meson loops should reduce the strength of spontaneous 
chiral symmetry breaking and, thus, reduce the values 
of | (q q) | and f n . Indeed, this is what has been found in 
Refs. |28l . I29I ] for the local NJL model. However, as seen 
in Figs. [4] and [H in the nonlocal model we find exactly the 
opposite behavior for all sets of model parameters. Obvi- 
ously, this difference between the local and the nonlocal 
model calls for some clarification. 

As pointed out before, in the nonlocal model, after 
introducing the Gaussian form factor, all diagrams at 
any order are automatically finite. This is different 
from the local NJL model where, because of its non- 
renormalizability, it is necessary to introduce indepen- 



dent cutoff parameters for pure quark loops and for 
meson-quark loops, respectively. In [2^, [20] a Pauli- 
Villars rcgularization has been used for quark loops and 
a three-momentum cutoff Am for meson-quark loops. In 
order to study the transition from the nonlocal model to 
the local one let us construct a nonlocal model with three 
parameters 

1. parameter of nonlocality A 

2. parameter of quark loop regularization A q 

3. parameter of meson loop regularization Am 
The local model corresponds to the limit 

A -j. co , A g = Ag* , A M - A M , (25) 

while the nonlocal model without rcgularization can be 
obtained by setting 

A = A fit , A q -> 00 , A M -> 00 . (26) 

For definiteness, let us compare the local model [HI, [29|] 
with the nonlocal one from [l3j] with parameterizations 
fixed in the mean-field approximation. Note that for the 
given parameterizations, the MF quark condensates in 
the local and the nonlocal models agree within less than 
0.5 %. 

The next step is to consider the 1/N C corrections and 
to investigate the role of the mesonic three-momentum 
cut-off Am- For this purpose it is very instructive to 
study the ratio of the full quark condensate to the MF 
contribution (qq) / (qq) MF ■ In Fig. [8]wc compare the Am 
dependence of this ratio for the local NJL model as given 
in Ref. [28j] (dash-dotted line) to that of the nonlocal 
model (bold solid line) and its local limit (dashed line). 4 
It is very interesting that in the region below ~ 1.5 GeV 
these models predict a negative sign for the 1/N C correc- 
tion whereas for large mesonic cut-off the sign is positive. 
However, in the nonlocal model the absolute value of the 
correction saturates for Am larger than ~ 2.5 GeV, which 
is well above actual parameterizations for A q and Am in 
Ref. [29||. In fact, in Ref. (28| it was found that already 
for Am ~ 1250 MeV, the pion propagator gets unphysical 
poles, which was interpreted as a breakdown of the ex- 
pansion. In our present model no such unphysical effects 
are observed. 

In order to study the dependence of the sign of the 
1/N C correction to the quark condensate on the form fac- 
tor we consider the Lorentzian-type form factor f[v%) = 
1/(1 + (p 2 E /A 2 ) n ) with n=2,5 or 10, see also (MEl- We 



4 Note that in Ref. [28l the total quark condensate was calculated 
including the perturbative part, whereas in the nonlocal model 
and its local limit we have only considered the nonperturbative 
part. However, for the local model the difference for nonstrange 
condensate is small. 
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FIG. 8: (Color online) The ratio (qq)/(qq) MF as a function 
of the meson cutoff Am- The full nonlocal result is shown by 
the black solid line. The local result (blue dash-dotted line) 
is taken from Fig. 3a of Ref. The "local" check (red 

dashed line) is the local limit of the nonlocal calculations, see 
Eq. (US]). 
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FIG. 9: (Color online) The ratio (qq)/{qq) MF as a function 
of the quark condensate for different form factors: Gaussian 
(G) and Lorentzian-type (L2, L5, L10) with n=2, 5, 10. For 
details, see text. 



found that the sign of the 1/N C correction is positive for 
all possible parametcrizations. This is shown in Fig. [9] 

A related question concerns the 1/N C corrections to the 
quark dressing functions A(p 2 ) and B(p 2 ), defined by 



Sp 1 = A(p 2 )$ - B(p 2 



(27) 



or, cquivalcntly, to the wave-function rcnormalization 
function Z(p 2 ) = 1/A(p 2 ) and the mass function 
M(p 2 ) — B(p 2 ) / A(p 2 ). In mean-field approximation, we 
have A MF (p 2 ) = 1 and B MF (p 2 ) = m c + m d f 2 (p 2 ), sec 
Eqs. and ©. Including 1/N C corrections we find that 
the vector dressing function A(p 2 ) becomes enhanced, 
corresponding to a reduction of Z(p 2 ). However, the 
scalar dressing function B(p 2 ) increases even stronger, 
so that the mass function M(p 2 ) is enhanced as well. 
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FIG. 10: (Color online) Behavior of the quark mass func- 
tion (upper panel) and wave-function renormalization func- 
tion (lower panel) for set No. 4 (bold lines) and set No. 2 
(thin lines). The results of the full 1/N C corrected approach 
(black solid lines) and the mean field contribution (red dashed 
lines) are compared with quenched (open symbols) and un- 
quenched (full symbols) lattice data }18j |. extrapolated to the 
chiral limit. 



In Fig. [TU] we present a comparison of the behavior 
of the mass function (upper panel) and wave-function 
rcnormalization function (lower panel) in the present ap- 
proach for two parameter sets (No. 2 and No. 4) with 
lattice data for this quantity [HI extrapolated to the chi- 
ral limit case. One can see that the gross behavior of 
the momentum dependence of the quark mass function 
as measured in lattice QCD is reflected by the nonlo- 
cal model. More in detail, the 1/N C corrected result for 
the weaker coupling No. 2 (thin black solid line) matches 
the quark mass at zero-momentum, but stays above the 
lattice data for dynamical quarks at finite momentum 
transfer. The stronger coupling set No. 4 (bold black 
solid line), however, overestimates the quark mass func- 
tion for low momenta p 2 < 0.5 GeV 2 by about 25% and 
stays below the lattice data by about the same amount in 
the intermediate momentum range 0.5 < p 2 [GeV 2 ] < 2. 

As one can see in the lower panel of Fig. \TU[ the agree- 
ment of our model results for the wave-function renor- 
malization function Z(p) with the lattice data is rather 
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is strictly positive. Hence, the pion gives a positive con- 
tribution to B N "^(p 2 E ), while the sigma contribution 
is negative because of the extra minus sign. Since the 
pions are more important, both, because of their lower 
mass and because of the degeneracy factor of 3, the total 
contribution is positive, as we have seen in Fig. 111! 

Neglecting the form factors and the momentum de- 
pendencies of the quark masses and coupling constants, 
Eq. (|28|) has a natural interpretation in terms of a simple 
model where the quarks are coupled to pions and sigma 
mesons by a Yukawa interaction, 

^Yukawa = —q{x) [g a <j{x) + g^ijsr ■ 7?(x)] q(x) . (30) 



FIG. 11: (Color online) l/iV c -correction to the scalar quark 
dressing function B(p 2 ) at p 2 = as a function of the dimen- 
sionless coupling GA 2 . Also shown are the separate contribu- 
tions of the diagrams (a) and (b) of Fig. O 



poor. This is not unexpected because Dyson- Schwinger 
QCD calculations find a 30% reduction of Z(p) at low mo- 
menta already in quenched approximation, i.e., without 
including meson effects [43|, whereas the present model 
predicts Z = 1 at mean field. A nontrivial mean-field 
Zip) can be obtained in the nonlocal PNJL model by 
adopting a suitable separable form of the interaction with 
vector currents with derivative coupling (rank-2 interac- 
tion) , see pjl where also a comparison with lat- 
tice QCD data can be found. In the present work which 
uses a separable ansatz in the scalar-pseudoscalar channel 
(rank-one interaction) , a nontrivial wave- function renor- 
malization arises only from the 1/N C corrections. Since 
these contributions generally do not exceed the 10% level 
when compared to the mean-field, we cannot expect a 
quantitative description of Z(p) at low momenta. 

In Fig. [Tl] we show the \/N c correction to -B(O) as 
a function of the dimensionless coupling GA 2 . Besides 
the total correction (black solid line), which is almost 
independent of the parameter choice, the contributions 
of the diagrams (a) and (b) of Fig. [3] are also indicated 
separately. As one can see, the dominant contribution 
comes from diagram (b), which is always positive. 

Indeed, in Euclidean space the contribution of diagram 
(b) to the i?-function takes the form 



d A q E P{ql)m{q\) 
(2^)4 q\ +77^2) 



2 

1e 



Ml ql + Ml 
where we have defined q+ = qE + Pe and 



M 2 M 



G- 



n M (g|) 



(28) 



(29) 



The latter is a straightforward extension of Eqs. (J7J) and 
(ITU|) to (Euclidean) off-shell momenta. In particular, g\j 



From this point of view, the observed enhancement of the 
£?-function through pion loops should be the expected 
result. This seems to be in conflict with the Dyson- 
Schwinger analysis of Ref. [43j where pion loops give a 
negative contribution to the _B-function in QCD. How- 
ever, in that approach the pions have been introduced 
somewhat differently, namely through corrections to the 
quark-gluon vertex. It is interesting that this leads to 
the opposite sign. Understanding this different behavior 
deserves further study. 

At first sight, our results also seem to contradict the 
analysis of Ref. [H[ in a local NJL model, where again 
it was found that the -B-function decreases when pion 
and sigma loops are included. In contrast to our present 
model, the l/A c -corrections have been included selfcon- 
sistently, i.e., the calculations take into account the back- 
reactions of the meson loops on the mean-field selfenergy. 
A closer inspection [47| reveals that this is the essen- 
tial difference. Whereas, in complete agreement with our 
expectations, the meson loops themselves give a posi- 
tive contribution to the i?-function, their back-reaction 
strongly reduces the Hartree contribution, so that the to- 
tal effect is a reduction of B(0). In turn, this suggests 
that the enhancement of B(0) we obtain in our model 
could be due to the non-selfconsistent treatment of the 
l/A c -corrcctions. On the other hand, it is a well-known 
problem of the fully sclfconsistcnt scheme that the in- 
ternal pions arc much too heavy, so that in Ref. [HI 
the (positive) contribution of the pion loops is underesti- 
mated. Hence, the true sign of the correction is the result 
of a delicate interplay between several different processes 
and certainly needs further investigations. 



III. FINITE TEMPERATURE 

In this section we extend the model to finite tempe- 
rature and then discuss the predictions for the pressure, 
the behavior of the quark condensate and critical tem- 
peratures. 
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A. Thermodynamic potential 

The model can easily be extended to finite temperature 
using a ^-derivable ansatz (see, e.g., Ref. [HI) supple- 
mented by the 1/N C expansion. The central quantity for 
the analysis is the thermodynamic potential per volume 



n = iTr MS- 1 ) + iTr(SS) + *(S) + £/($, $) - fi , 



where S and £ = (S' c ) _1 — S _1 are the full propagator 
and the quark selfenergy, respectively, and Tr denotes 
the trace over all degrees of freedom, internal ones and 
4-momenta. At nonzero temperature, we also take into 
account the Polyakov-loop dynamics which no longer de- 
couples from the quark sector. To that end a constant 
temporal background gauge field <j> = (A4) = (iAo) is 
minimally coupled to the quarks and a Polyakov loop 
potential is added in Eq. (HU). Here $ = 

w- Tr c e % $l T denotes the Polyakov loop expectation value 
and $ its conjugate. In order to avoid confusion with the 
potential of the ^-derivable scheme ( "$ functional" ) , we 
denote the latter as ty. In the exact case, "J contains all 
two-particle irreducible diagrams. Finally, wc have intro- 
duced a subtractive renormalization constant f2 , which 
is chosen such that the vacuum (T = 0) has vanishing 
pressure. 

The thermodynamic equilibrium corresponds to the 
(global) minimum of the thermodynamic potential with 
respect to the full quark propagator and to the Polyakov 
loop, so that the following necessary conditions (gap 
equations) must be fulfilled 



an 

dS 



0, 



an 
9$ 







on 



o 



(32) 



We work in Polyakov gauge where the background gauge 
field is dia gon al in color space, i.e., <fi = 4*3X3 + ^sAs- 
Following [1J|, we require $ = <I> to be real with real 
03, </>8- As a consequence 4>s = and we are left with one 
variable (/) 3 . The second and the third equation thus re- 
duce to dSl/84>3 = 0. Moreover, the first equation implies 
that £ = iS^/SS. Diagrammatically this means that the 
quark selfenergy can be obtained from \P by cutting one 
quark line at all possible places. 

Approximations can be introduced by truncating \1/ 
at a certain order. The mean-field results, correspond- 
ing to the leading order in l/N c , are obtained from the 
"glasses" diagram, displayed in Fig. [T2] Solid lines rep- 
resent dressed quark propagators. In the nonlocal model 




FIG. 12: Glasses diagram for the W potential in the $- 
derivable scheme. 



+ 
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FIG. 13: Ring sum in the ^-potential, see Eq. (|35|) . 



(31) the "glasses" potential takes the form 



= E % [-iv(r M *s)f, 



glasses 



(33) 



and the thermodynamic potential reads 



2G 



where k"i 



i=0, ± 



(oj l n ) 2 + k 2 and the notation J 



TJ2J d 3 k/(2n) 3 has been introduced. Note that due 

n 

to the coupling to the Polyakov loop the fcrmionic Mat- 
subara frequencies uj n — (2n + 1)ttT are partially shifted: 
uj^ = u n ± 4>3 , u>® = u) n . Apart from this shift, the 
modification of the quark propagator only depends on 
the dynamical mass md- Therefore at mean field the gap 
equations, Eq. (f32|) . take the form 



MF 



dm d 



0. 



on 



MF 



0. 



(34) 



The next-to-lcading order contribution to the ']>- 
potential is given by the "ring sum" , 



E ^THn[l-GII M ], 

M = 7T.<T 



see Fig. [13] Here II M denotes the quark-antiquark polar- 
ization functions constructed with the full quark propa- 
gators and d,M is the mesonic degeneracy factor. 

At this level, the problem arises that in a fully self- 
consistent treatment the iteration of diagrams in the 
gap equation leads to contributions of arbitrary orders 
in l/N c . As a consequence, different approximation 
schemes can be defined. In the present paper, we use a 
"strict 1 /N c expansion" , where all contributions beyond 
the next-to-leading order are discarded. In the absence 
of the background gluon field this scheme is straightfor- 
wardly implemented by first solving the mean-field gap 
equation for nid (first equation of Eq. ([34])) and then 
evaluating the ring sum using the mean-field propaga- 
tors. Thus, one gets for the thermodynamic potential 



O 



MF 



(36) 



with the 1/N C correction 

qNc = [ ln[l-GH M (p> m )] (37) 



M= 



10 



where the polarization functions IT M are evaluated with 
mean-field propagators and summation is over bosonic 
Matsubara frequencies. 

Including the gluon background, we suggest to treat 
the Polyakov-loop potential as effectively N c indepen- 
dent. 5 A strict 1/N C expansion of the thermodynamic 
potential then corresponds to evaluate Eq. for the 
simultaneous solutions of the gap equations 



MF 



dil 

dm d d(f>3 



0. 



(38) 



Note that 03 is determined by minimizing the total ther- 
modynamic potential, whereas m,d is obtained from the 
mean-field part only. Nevertheless, since f2 MF also de- 
pends on 03, the value of md is changed as well compared 
to the mean-field calculation, due to the modified value 
of 3 . 

We also note that the scheme outlined above slightly 
differs from the scheme in our previous paper [l3j where 
both, nid and 03, have been fixed at mean-field level, 
before inserting them into Sl Na . In the numerical calcu- 
lations, however, this difference turned out to be small. 



B. Finite temperature results 

We now discuss our numerical results at finite tempe- 
rature. In the quark sector we take the parameters of set 
No. 4, see Table [J in Sec. HTCl For the Polyakov loop 
potential U($, $) we adopt the logarithmic parametriza- 
tion of Ref. [1JI, which has been fitted to the quenched 
lattice data of Ref. [5l| . The only exception concerns the 
parameter To, which corresponds to the transition tem- 
perature in pure gauge. While in Ref. fioj the empirical 
value Tn = 270 MeV was used, it has been argued in 
Ref. [52| that To depends on the number of active fla- 
vors and a value of 208 MeV was suggested for Nf = 2. 
Therefore our calculations will be performed using both 
values of To • 

The model predictions for the pressure are displayed in 
Fig. [H The upper panel corresponds to T = 270 MeV, 
the lower one to T = 208 MeV. Together with the full 
result (solid lines) we show the partial contributions to 
the pressure from the mean-field approximation (dash- 
dotted lines), P MF = -n MF , and from the 1/N C correc- 
tions (dashed lines), P Nc = — S! Nc . In agreement with 
earlier results in PNJL-like models we find that at low 



5 In principle the [/(<t>, <t>) is proportional to the number of gluons, 
N% — 1 |49l. l50ll . Its leading contribution to the thermodynamic 
potential is therefore of the order O(N^), while the quarks only 
contribute at the order 0(N C ) and corrections are of the order 
0(N®) for both, quarks and gluons. However, since in practice 
the detailed form of U is not based on a 1/N C expansion, but 
rather a phenomenological parameterization fitted to quenched 
lattice data, we believe that it is more appropriate to treat it as 
N c independent in the present context. 
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FIG. 14: (Color online) Temperature dependence of the scaled 
pressure P/T 4 for T = 270 MeV (upper panel) and for To = 
208 MeV (lower panel): mean field contribution (blue dash- 
dotted lines), 1/N C correction (red dashed lines), mean field 
+ 1/N C contributions (black solid lines). The green dotted 
lines correspond to the NLO chiral perturbation theory (xPT) 
result, Eq. (O [H. 



temperatures the mean-field contribution, corresponding 
to thermally excited quarks, is strongly suppressed by 
the Polyakov loop. In this regime the thermodynamics is 
governed by the ring sum, which is dominated by pionic 
degrees of freedom as the lightest particles in the mass 
spectrum. Therefore, it is instructive to compare our re- 
sult with the predictions of chiral perturbation theory 
(xPT). xPT describes the low-energy structure of dif- 
ferent amplitudes in terms of an expansion in powers of 
energies, momenta and current quark masses. The finite 
temperature result for the pressure is given by [53| 



P. 



X PT 



N]-l 



9o 



1 Ml 9l ) + 0(P S ) , 



fjo 
9i 



-IT 



N f 2F 2 
d 3 p 



(2tt) 3 



d 3 p 



In (l - e- E -' T 
1 



(2tt) 3 E„ (e s -/ T - 1 



(39) 



where F is the weak pion decay constant in the chiral 
limit, En = \/p 2 + A/2 is the pion energy, and 0(p 8 ) 
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refers to the chiral counting scheme, where M v and T 
count as quantities of order p. The term proportional to 
go just represents the free relativistic pion gas pressure, 
while the g\ term is caused by interactions and leads to 
a small reduction of the pressure. 6 The omitted terms of 
order 0(p 8 ) are also due to interactions. 

Recalling that the pion decay constant is of the order 
we see that the go and g\ terms are of the order N® 
and 1/N C , respectively. Comparing this with our model, 
where the mean-field and ring-sum contributions to the 
pressure are of the order N c and N®, respectively, we 
conclude that our model calculations should only be con- 
sistent with the lowest-order X-PT result, i.e., with the 
go term. Nevertheless, since the g^-term is small (it van- 
ishes in the chiral limit), we find excellent agreement even 
when this term is included, see Fig. Q31 Moreover, the 
low-temperature behavior of our model predictions is al- 
most insensitive to the particular functional dependence 
of the form factor and different parameterizations. Of 
course, our results start to deviate from the xPT predic- 
tions when we approach the chiral phase transition. Near 
the critical temperature the a meson gives an additional 
visible contribution whereas already for T > 1.5 T c the 
mesonic contributions are negligible and the quark-gluon 
mean-field dominates the pressure. 

Comparing the upper and lower panels of Fig. [H] we 
observe that the lowering of To also leads to a lowering 
of the transition temperature. However, apart from this 
trivial effect, the results remain qualitatively unchanged. 

In Fig. [T5] we show the temperature dependence of 
the quark condensate (qq)r (black solid line) and of the 
Polyakov loop expectation value (red dashed line) in our 
model beyond mean field. For comparison we also show 
the quark condensate in mean-field approximation (blue 
dash-dotted line). Again, the results for T = 270 MeV 
and 208 MeV are displayed in the upper and in the lower 
panel, respectively. As already seen for the pressure, the 
1/N C corrections mainly affect the behavior of the chiral 
condensate below and around the critical region. In par- 
ticular at low temperatures the reduction of the chiral 
condensate is almost entirely driven by the pion dynam- 
ics, which is missing in the mean-field approximation. In 
this context, it is again instructive to compare the model 
results with the X-PT prediction. The latter is given by 



(qq) i 



i 



9i 



Nf F 2 



o( P 4 



(41) 



showing that the leading temperature effect is of the 
order 0(1/N C ). Our model results are completely con- 
sistent with this: Whereas the mean-field condensate 



In the chiral limit go and g\ are given by 

7T 2 T 4 T 2 

1 91 = — • 

45 12 



90 



(40) 



(0(1/ N®)) stays practically constant at low T, we find 
very good agreement with the xPT predictions (green 
dotted line) after including the 1/N C corrections. We 
should note that the chiral expansion scheme is very sim- 
ilar but not exactly equivalent to the 1/N C counting: In 
\PT the leading-order correction to the chiral conden- 
sate depends on the weak pion decay constant in the chi- 
ral limit F, whereas in the 1/N C expansion it depends on 
the week pion decay constant for massive pions in mean- 
field approximation /^ F . Although formally of higher 
order in p or N c , this could lead to quantitative differ- 
ences, even at low temperatures. However, for the chosen 
parameters the difference between F and f^f happens 
to be very small. 

In agreement with other works [29|, |33|, HH , the addi- 
tional reduction of the chiral condensate by the 1/N C - 
correction terms generally also leads to a lowering of the 
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FIG. 15: (Color online) Temperature dependence of the quark 
condensate normalized to its vacuum value A = (qq)r / (qq) 
(black solid lines) and the Polyakov loop (red dashed lines) 
in the nonlocal PNJL model beyond mean field. Furthermore 
shown are the mean-field result for the normalized quark con- 
densate A MF (blue dash-dotted lines), the lowest-order chiral 
perturbation theory (\PT) result (green dotted lines), and a 
naive polynomial interpolation in the unstable regions of the 
1/N C expansion (thin black solid lines). The model calcula- 
tions have been performed with To = 270 MeV (upper panel) 
and with To = 208 MeV (lower panel). 
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chiral phase-transition temperature as compared to the 
mean-field result. This effect is most clearly seen in the 
upper panel of Fig. Q~5] corresponding to To = 270 MeV. 
Unfortunately, because of the perturbative nature of the 
strict l/A^-expansion scheme, our model cannot be ap- 
plied to study the phase transition itself. In the figure, 
this is obvious from the existence of an unstable region 
("wiggle") in that region. The origin of this wiggle can 
be attributed to the momentum independent 1/N C cor- 
rection to the quark selfenergy (Fig. In the vicinity 
of the phase transition, this diagram is dramatically en- 
hanced due to the restoration of chiral symmetry and 
the corresponding lowering of the intermediate er-meson 
mass. (In fact, in the chiral limit this contribution would 
go to minus infinity at T c .) As a rough estimate we de- 
fine the unstable region as the regime where the relative 
correction to the quark condensate is larger than 1/N C . 
In the present example this corresponds to temperatures 
between 183 and 223 MeV. Since our results cannot be 
trusted in this area, we suggest to use a simple poly- 
nomial interpolation between the stable regions at lower 
and higher temperatures. The resulting temperature de- 
pendence is displayed in Fig. [15] by the thin black solid 
line. 

The situation is somewhat different for To = 208 MeV 
(lower panel). Although the 1/N C corrections are again 
essential at low temperatures (in agreement with the 
xPT results), their effect on the chiral critical tempera- 
ture is almost negligible. This is due to the fact that for 
this lower value of To the mean-field transition takes place 
earlier. Related to this, the mean-field chiral condensate 
drops more steeply 7 . The "unstable region" where the 
relative correction to the quark condensate is larger than 
1 /N c is quite narrow and the wiggle near the chiral phase 
transition, which we found for To = 270 MeV, is not 
present for T = 208 MeV. 

A systematic overview about the To-dependence of 
the pseudocritical temperatures for deconfinement (red 
dashed line) and chiral restoration (solid black line) is 
given in Fig. [T75J As definitions of these temperatures we 
have chosen the maxima of the temperature derivatives 
of the corresponding order parameter, i.e., Polyakov-loop 
expectation value and chiral condensate, respectively. In 
the unstable regime (indicated by the shaded area) the 
latter was obtained from the polynomial interpolation. 
The mean-field result for the chiral-restoration tempera- 
ture is also shown (blue dashed-dotted line). 

In this figure, we can roughly distinguish three regimes: 
At large values of To already the mean-field T c is lower 
than Td, and T c is even lowered further by including 1/N C 
corrections. Then, in some intermediate regime, Td and 
T c agree at mean field but there is still a visible reduction 
of the latter by the 1/N C corrections. Finally, for Tq < 



7 For To < 150 MeV, the phase transition even becomes first order. 
A similar observation was recently made in Ref. [lBl . 
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FIG. 16: (Color online) Critical temperatures versus the pa- 
rameter To of the Polyakov-loop potential: deconfinement 
(red dashed line), chiral restoration in mean-field (blue dash- 
dotted line), chiral restoration with mesonic fluctuations 
(black solid line). The region with | (g<7)T7(g<7)T F > 1/N C 
is indicated by the shaded area. 



220 MeV all three lines practically coincide. 8 

To understand this behavior we should recall that 
without coupling between quark and Polyakov-loop sec- 
tor, To corresponds to the (first-order) deconfincmcnt- 
transition temperature, while the chiral cross-over in 
the nonlocal model takes place at a very low tempera- 
ture, T C MF = 116 MeV. From this point of view it is 
quite remarkable that even for To = 270 MeV the cou- 
pling reduces the difference between T C MF and Td to less 
than 10 MeV. Nevertheless, the synchronization obvi- 
ously works better when To is reduced and eventually 
becomes perfect for T < 220 MeV. Related to this, the 
transitions get sharper with decreasing To and the rel- 
ative effect of the 1/N C corrections on T c is reduced as 
well, as we have seen above. 

At this point we should note that the Polyakov-loop 
dynamics and, hence, the deconfinement transition are 
quite insensitive to the 1/N C corrections. This is most 
likely an artifact of the model where the Polyakov-loop 
is affected by the meson back-coupling only in a very 
indirect way (see the discussion after Eq. ([38])). In prin- 
ciple, mesonic fluctuations should also directly influence 
the Polyakov-loop potential. However, it is not clear how 
to include this in PNJL-like models, where the potential 
is obtained from a phenomenological fit of quenched lat- 
tice data where even the effects of dynamical quarks arc 
neglected (except for changes in To). In this context it 



According to Table [I] the 1/N C corrections can even enhance the 
value of T c by 1 or 2 MeV at To = 208 MeV. However, since cross- 
over temperatures are not defined in a unique way, and even rely 
on the polynomial interpolations in the present case, this should 
not be over-interpreted. If we had defined T c as the temperature 
where the condensate reaches half of its vacuum value, the 1/N C 
corrections would always lead to a reduction of T c . 



13 



1.0 



e 0.4 

<f 

0.2 
0.0 



■ 


A 




A 

$ 


*\\ 

M 




A ltS , 


hisq a / 




^reni 


hisq *y 






stout T, . 


if 


1 x 1 L ien , 


stout j& 
■ 4* ^ 





50 100 150 200 250 300 350 



T [MeV] 



is 0.8 and the ratio of the nonstrangc current quark mass 
to the strange one is about 1/20. On the other hand the 
main role for low temperature changes of the quark con- 
densate is played by pions - an effect which is correctly 
accounted for in our model. 

Indeed, we find that A; iS is rather well reproduced by 
our model calculations for A. For the Polyakov loop, 
on the other hand, we find good agreement only at the 
onset, whereas around the dcconfinement temperature, 
T w 170 MeV, the model prediction rises much steeper 
than the lattice data. This seems to indicate that there is 
missing physics, which is not captured by just rescaling 
the temperature Tq of the pure-glue potential. 



FIG. 17: (Color online) Temperature dependence of the quark 
condensate normalized to its vacuum value A = (qq)T/{qq) 
(black solid line) and the Polyakov loop (red dashed line) in 
the nonlocal PNJL model with To = 208 MeV beyond mean 
field in comparison with lattice QCD data for the subtracted 
quark condensate A( i3 and for the renormalized Polyakov 
loop I/rcn- Lattice points are taken from [l9(] (hisq action - 
HotQCD collaboration) and [STJ] (stout action - Wuppertal- 
Budapest collaboration). 



would certainly be interesting to compare the behavior 
of the Polyakov loop with the so-called dressed Polyakov 
loop, related to the dual quark condensat e loll IHrjj . The 
latter can be calculated within PNJL-like [57| (and even 
NJL models (HI) and it should be possible to study l/N c 
corrections as well. However, this is beyond the scope of 
the present paper. 

In Fig. [17] wc present the comparison of our predic- 
tion for the temperature dependence of quark condensate 
and Polyakov loop with recent lattice data [l9l . [20j . The 
model calculations are performed with To = 208 MeV, 
which leads to the more realistic transition temperatures. 
One should mention that the lattice calculations [l9|, [2(3] 
are performed for the case of two light quarks and one 
strange, Nf = 2+1. For the investigation of the quark 
condensate they use a combination of nonstrange(Z) and 
strange(s) quark condensates, called subtracted quark 
condensate, 



(42) 



In the present paper we consider a model with only two 
light active flavors. Nevertheless, it makes sense to per- 
form the comparison of our calculation for the tempera- 
ture dependence of the quark condensate normalized to 
its vacuum value A = (qq) t /(?<?) with the lattice calcu- 
lations of the subtracted quark condensate. The reason 
is that the nonpertubative part of the quark condensate 
in Eq. ([6]) is finite and the strange quark contribution 
in Eq. (|42p is suppressed in the nonlocal model by two 
factors: (i) the ratio of the nonpertubative part of the 
strange quark condensate to that of the nonstrange one 



IV. SUMMARY 

In the present work, we have extended the nonlocal chi- 
ral quark model, coupled to the Polyakov loop, beyond 
the mean field approximation using a strict 1 /N c expan- 
sion scheme. In vacuum, it is found that 1/N C corrections 
lead to an increase of the quark condensate, which is op- 
posite to results obtained in the NJL model with local 
interactions. However, the result of the local model is 
strongly dependent on the mesonic cut-off and for large 
values of the cut-off (Am > 1-5 GeV) the sign of the 
correction is positive in both models. 

The parameters of the nonlocal model have been re- 
fitted in vacuum to reproduce the physical values of the 
pion mass and the weak pion decay constant after in- 
cluding the 1/N C corrections. In agreement with general 
expectations, we find that the mean field gives the domi- 
nant contribution to the pion properties, while the max- 
imal size of the 1/N C corrections to the pion mass and 
weak pion decay constant amounts to 15 and 20 MeV, 
respectively. 

At finite temperature, the 1/N C corrections lead to a 
reduction of the chiral condensate when compared to the 
mean- field result. Typically, this also leads to a reduc- 
tion of the chiral phase transition temperature. However, 
for lower values of the parameter To in the Polyakov-loop 
potential, the mean-field transition becomes steeper and, 
thus, the effect of the 1/N C corrections on T c becomes 
smaller. Eventually, for To < 220 MeV, T c remains prac- 
tically unaffected by the 1/N C corrections. For low tem- 
peratures, T < 100 MeV, our result for the quark conden- 
sate practically coincides with that of X-PT whereas the 
high-temperature region is well controlled by the mean- 
field approximation. 

In the chiral limit our expansion scheme, which treats 
the l/N c corrections perturbativcly, breaks down in the 
vicinity of the chiral phase transition. For the real 
case of nonzero current quark masses a similar situa- 
tion takes place for large values of the To parameters 
of the Polyakov loop potential (e.g., To = 270 MeV in 
the pure gluodynamics case). Alternatively, one could 
include the mesonic fluctuations non-perturbatively, as 
done in Ref. [3J for the local NJL model. In that anal- 
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ysis it was found as well that mesonic fluctuations lead 
to a decrease of the critical temperature, but the un- 
stable region around the phase transition is absent. On 
the other hand, the low-temperature behavior of xPT is 
not reproduced in this scheme. It seems plausible that 
the inclusion of higher-order corrections could make both 
expansion schemes converge to one another, i.e., in the 
strict 1 /N c expansion the unstable region will be reduced 
(except close to a real phase transition) while in the self- 
consistent scheme the low-temperature behavior will get 
closer to the xPT predictions. 

Concerning the pressure, we confirm our previous re- 
sult (l3j that the pionic contribution dominates in the low 
temperature region. In this regime, the pressure is quite 
insensitive to the details of the interaction and agrees 
almost exactly with that of an ideal pion gas. At tem- 
peratures T > T c , on the other hand, the mesonic con- 
tributions die out at T ~ 1.5 T c . 

As a next step, we plan to study nonzero chemical po- 
tentials, which may require a nontrivial extrapolation of 
the Polyakov-loop potential into this regime [52j - The 
treatment of the entire T-\i plane of the QCD phase dia- 
gram finally requires the inclusion of baryonic degrees of 
freedom as well. 
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where the factor 4 in II^3 Q (p 2 ) is the degeneracy factor. 

The calculation of the weak pion decay constant is 
more complicated. Namely, part of the diagrams can be 
obtained from the pion propagator times — ig^ F and sub- 
stitution of the outgoing pion vertex by the vertex with 
an external current. Furthermore, there are additional 
nonlocal diagrams, shown in Fig. 1191 

The additional 1/N C corrections of type 1 are 
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Appendix A: l/N c corrections to the pion 
propagator and the weak decay constant 

The 1/N C corrections to the meson propagator consist 
of the diagrams shown in Fig. [T8J For convenience we 
divide them into three types. For the pion the corre- 
sponding expressions take the form 

d A k 



n 



Nc (v 
x Tr 



n 



N c / 2\ 

7r,lb+2hlP ) 



(2^) 4 

-p7T o yiN c Q p7T Q 

1 fc + ,fc„' : 'fc- Zj fc_ D fc- i k-,k+ D k+ 

d 4 k 



(27TY 



x Tr 



the 1 /N c corrections of type 2 take the form 
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and should be summed over M = n, a. 
The additional type-3 corrections are 
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la lb 2a 2b 2c 3a 

FIG. 18: Diagrams for the calculation of 1/N C corrections to the meson propagator. 
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3b 3c 
FIG. 19: Additional nonlocal diagrams for the calculation of 1/N C corrections to the weak pion decay. 
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